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The casing of large electrical generators can be deformed slightly by the rotor’s mag- 
netic field. The sound emission produced by these periodic deformations, which could 
possibly exceed guaranteed noise emission limits, is analysed analytically and numeri- 
cally. From the deformation of the casing, the normal velocity of the generator’s surface 
is computed. Taking into account the corresponding symmetry, an analytical solution 
for the acoustic pressure outside the generator is found in terms of the Hankel func- 
tion of second order. The normal velocity of the generator surface provides the required 
boundary condition for the acoustic pressure and determines the magnitude of pressure 
oscillations. For the numerical simulation, the nonlinear 2D Euler equations are formu- 
lated in a perturbation form for low Mach number Computational Aeroacoustics (CAA). 

The spatial derivatives are discretized by the classical sixth-order central interior scheme 
and a third-order boundary scheme. Spurious high frequency oscillations are damped by 
a characteristic-based artificial compression method (ACM) filter. The time derivatives 
are approximated by the classical 4th-order Runge-Kutta method. The numerical results 
are In excellent agreement with the analytical solution. 


Introduction 

For large air-cooled electrical generators, periodic 
deformations of the generator casing are induced by 
the rotor’s magnetic field (Fig. 1). Finite element 
based structural analysis suggests deformation ampli- 
tudes relative to the generator radius less than 10“ 5 m , 
cf. Ahrens 1 . Albeit these deformations are quite 
small, they are a source of sound. They could pos- 
sibly exceed guaranteed noise emission limits. In this 
note we present analytical and numerical solutions of 
the sound pressure level resulting from rotor induced 
casing deformations with given amplitude. 
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Analytical Solution 

To simplify the problem, we consider a generator 
of infinite axial extent. The problem then reduces to 
two dimensions (in the x : y or the r, <f> plane, respec- 
tively). The deformation of the casing induced by the 
rotor's magnetic field provides a (velocity) boundary 
condition for the acoustic field. Note that oniy the ve- 
locity component normal to the surface of the casing is 
relevant for the acoustic field. Therefore the acoustic 
pressure distribution produced by the casing deforma- 
tions is equivalent to that of a rotating ellipse (Fig. 
2 ). 

As a first step we determine the velocity of the sur- 
face of the ellipse resulting from its rotation. Without 
loss of generality, an ellipse oriented as shown in Fig. 
2 at time t = 0 is considered. The points x s on the 
surface of the ellipse may be described by 

( a cost \ 

6 sin t )' 

where a, b are the semi-major and semi-minor axes, 
respectively, and the parameter r runs from 0 to 27r. 
We assume that the casing deformations are small in 
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Fig. X Schematic of the casing deformation induced by the magnetic field of the rotor. The deformation 
of the casing follows the rotor’s movement. 


amplitude, and write 

CL = .R(l d" *0) 

b = R{ l-€), 

with e <£C 1. The following approximations are then 
accurate to second order in e: 

tan 0 % (1 — 2c) tan r, 

r 6 4- e sin 20, 

where 0 is the azimuth angle, as indicated in Fig. 2. 

Using Taylor expansions and trigonometric identi- 
ties, it is straightforward to show that at time t — 0 
the surface of the ellipse is described by 

x s (^,t = O) = R(l+ecos20)( )+C?(f 2 * ). (1) 

At later times £, the surface vector x. s at angular po- 
sition 0 is equal to 

x 5 (0, t) = 7l(Qt) x s (0 — Qt, t = 0), (2) 

where 

IJ/m _ f COS0 - S ' m0 \ 

^-{sind cos 9 ) 

is the transformation matrix that rotates a vector by 
an angle 9 about the origin in the (x,y) plane. Equa- 
tion (2) expresses the fact that the ellipse rotates about 
its center by an angle equal to Vtt during time t. Com- 
bining (1) and (2), we obtain 

x,0M) = /?[l+fcos2(<?-fi<)] ( ) + 0 (c 2 ). 

(3) 

Differentiating with respect to time yields the velocity 

V.(*,t) = 2«flnsin2(*-ni)f )+C?(f 2 ). (4) 

Note that as we consider the movement of points on 
the ellipse’s surface along a given direction ( we keep 0 


fixed), the velocity vector v s {<p,t) is aligned with the 
radial direction. As the ellipse’s surface normal n = 
(cos 0; sin 0) 4- O(e) also points approximately in this 
direction, we obtain for the surface normal velocity 

v s ■ n(0 , t) zz Vf, r (<p, t ) = 2 eR H sin 2(0 — Qt) + C>(e 2 ) 

for the radial velocity component v s>r . Introducing 
the more convenient complex exponentials instead of 
the sine function, where it is understood that physical 
quantities are the real part of complex numbers, vre 
obtain 

i; 3t r (<M) = -if/? u>e 2i *e- iut + 0{ f 2 ), (5) 

where oj = 2Q has been introduced. We find that 
according to (3) and (5), each surface element of the 
generator oscillates harmonically with amplitude eR 
and a frequency u } which is equal to twice the angular 
frequency Q of the rotor. 

Because the displacement of the generator surface 
is small, it is appropriate to solve for the distribu- 
tion of pressure disturbances outside a circular disk 



Fig. 2 Schematic and relevant notation of the 
equivalent problem of a rotating ellipse. 
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of radius R and enforce the boundary condition (5) 
at r = R. The acoustic fluctuations of pressure p' 
obey the Helmholtz equation, which reduces for the 
harmonic case under consideration to 

V 2 p' + k 2 p' = 0, (6) 


with the wave number k = oj/c (c is the speed of 
sound). After separation of variables 


one obtains 



(?) 


Due to the angular dependence given by (5), <$(<£) = 
exp{2i<£} and m = 2. The general solution for the 
radial function g(kr) is then (see Morse & Ingard, Sect. 
7.3 2 ) 

g{kr) = AH^ikr). (8) 

where = Jm + iN m is the Hankel function of 

mth order ( J m and N rn are Bessel and Neumann func- 
tions, respectively), and A is a constant which has to 
be determined from boundary conditions. At the (im- 
permeable) surface of the generator, 


is also described by (5). Hence the acoustic pres- 
sure distribution (10) produced by a rotating ellipse 
of small eccentricity is identical to that of a Kirchhoff 
vortex of equal size, eccentricity and rotational speed. 
Farfield and nearfield approximations of the pressure 
distribution (10) are also presented in Muller 3 , and 
it is pointed out that a farfield approximation derived 
using Green’s function yields identical results. 

However, for large electrical generators used in 
power generation, typical Helmholtz numbers Ti might 
range from 1 to 5. With Ti ~ 0(1), ”small-source” or 
" far- field” approximations are not applicable when the 
pressure distribution in the vicinity of the generator 
has to be determined. 

Numerical Solution 

We follow the approach outlined by Muller and 
Yee 16 for the Kirchhoff vortex sound. Before the ap- 
plication of the numerical scheme, we conditioned the 
nonlinear 2D Euler equations by reformulating the 
equations in a perturbation form for low Mach number 
Computational Aeroacoustics(CAA), cf. Sesterhenn et 
al. 4 . The unknowns are the changes of the conserva- 
tive variables with respect to their stagnation values. 
This formulation retains all nonlinearities and the con- 
servation form, while numerical cancellation errors are 
minimized compared to the original conservative laws. 
The Euler equations in this perturbation form can be 
written as 


and the linearized form of the radial momentum equa- 
tion 

dv! r _ _cy 

P dt Dt 

yields then with (5) the relations required to fix the 
constant A. We obtain the simple result 


A = 


peJ 2 R 
kH { 2 l) '(kR)’ 


( 9 ) 


where ” ,r denotes the derivative of a function with 
respect to its argument. 

The problem is non-dimensionalized conveniently 
with the generator radius R , the speed of sound c and 
the density p. The pressure fluctuations p' are non- 
dimensionalized with pc 2 . Introducing the Helmholtz 
number W = kR, we obtain finally 


P’(r, <t>) = exp{2^}. (10) 

where the tilde indicates non-dimensionalized quanti- 
ties. Numerical values for the derivatives ~ 

^2(2) exp{i72(2)} of the Hankel functions are given 
in Table V of Morse and Ingard 2 or can be computed 
with standard mathematical software. 

It has been shown by Muller 3 that the normal veloc- 
ity of the surface of an almost circular Kirchhoff vortex 


^ + v (pu)' = 0 , ( 11 ) 

+ V ■ (pu)'u' + Vp' = 0, (12) 
+ V ■ ((pH)'u' + (pH)on') = 0, (13) 

where 

P = p- Po, (pu)' = pu. ( pEY = pE- ( pE)o , 


(pu) 

+ P r 


p' = (7-l)[(^)'-i(pu)'.u'), 

(pHY = (pEY +p'. 


Here, p denotes the density, u the velocity, E the to- 
tal energy per unit mass, H the total enthalpy, and 
7 — 1.4 the ratio of specific heats for air at standard 
conditions. The ‘7 n and subscript U (T denote pertur- 
bation and stagnation variables, respectively. 

Although this perturbed form is identical to the 
original conservation laws, discretizing e.g. Vp leads 
to cancellation errors, whereas these errors are avoided 
when discretizing Vp'. In Cartesian coordinates, the 
perturbed 2D Euler equations can be expressed as 


01 y of ' 2 

dt Ox + dy 


(14) 
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where 


(*»)' \ 

(/**)'«' + p’ 

(pv)’u' 

(pH Yu' + (pH)ou' J 
\ 

‘ 2_ | wv+y 

V (Pif) V + (pF)ou' / 

Here, u* — u is the x-direction velocity and v l — v is 
the y-direction velocity. 

For the treatment of general geometries, a coor- 
dinate transformation (x(£, 7/),y(£,r;)) is used. The 
resulting transformed 2D Euler equations, transformed 
variables and fluxes are 


u 7 = 



F', = 


( (r»Y 

(pu)'v' 


di]' 0F\ 0F' 2 

j- H - 0 

dt d£ drj 


(15) 


U' = j-’u', 

F\ = + J-'SyF' 2 , 

F' 2 = J- x r, x F\ + J- l T, y F' 2l 


with the Jacobian determinant of the transformation 
J~ l = x^yrj — x T? y^ ; and the metric terms J~ l £ x — 

J — *£775 J Vx — ^ *]y X£, 

Kreiss and Scherer 5 , Strand 6 and Carpenter et 
al. 7 constructed high order difference operators Q for 
“e//dx” such that the summation-by-parts (SBP) prin- 
ciple holds, i.e. 


(u,Qv) h — uvvx - u 0 v 0 - {Qu,v) h , (16) 

where € R‘ v+1 , h is the uniform grid spacing and 
N 4- 1 is the number of grid points. The discrete scalar 
product and norm are defined by 

(u,v) u = hu T H n, ||«||fc = ( U,u) h , 


where FT is a symmetric positive definite (;Y + 1) x 
(N -1-1) matrix. 

We employ a Q operator, which is third-order ac- 
curate near the boundary and compatible with the 
standard sixth-order central difference operator in the 
interior. It was derived by Strand 6 and Q x , the x- 
direction Q operator for ^ is of the form 

{ £ djk v ki .7=0,..., 5, 

(Q^v) Jt j = 6, ...,N — 6, 

-\Yl 8 k=odv-j.k v N-k, j = N - 5, N , 

( 17 ) 

where (Q x v)j — h(go v >+3— 1 20^+2 + ~ 4 v j— 1 

— g^Uj_3) is the standard sixth-order central 
difference approximation of the first derivative. The 
forms of the 5x9 matrix D = (djk) and matrix H 


can be found in Strand 6 and Gerritsen and Olsson 8 . 
Here H is a diagonal matrix defining the norm of the 
Q operator. The global order of accuracy for (17) is 
four. Since (17) is based on a diagonal norm, its ap- 
plication to multi-dimensions is straightforward. The 
metric terms are discretized by the same difference op- 
erators as the flux derivatives in (15). In order to 
guarantee freestream preservation in 3D for high or- 
der spatial discretizations, the metric term treatment 
of Vinokur and Yee 11 is recommended. 

With the density, velocity and pressure non- 
dimensionalized with reference quantities po, Co and 
pocg, respectively, the Riemann invariants can be ex- 
pressed as p' - u n , p r — p* , u t and p' + u n . Here, 
u n is the normal velocity and u t is the tangential 
velocity. At the generator, the in-going Riemann in- 
variants are prescribed using the analytical solution 
(10), cf. Muller and Yee 16 . At the farfield, nonreflect- 
ing boundary conditions are implemented by setting 
the in-going Riemann invariants equal to zero cor- 
responding to stagnation conditions. The outgoing 
Riemann invariants are taken from the numerical solu- 
tion computed at the boundary. We use the so-called 
injection method to prescribe the in-going Riemann 
invariants, i.e. by imposing them explicitly which 
might destroy strict stability. Research is in progress 
to implement the in-going Riemann invariants by the 
penalty method called “simultaneous approximation 
term'’ (SAT) of Carpenter et al. 7 or the projection 
method of Olsson 1213 . Either approach yields strict 
stability for a Q operator satisfying the SBP principle 
(16) for linear problems. For nonlinear systems of con- 
servation laws, the in-going Riemann invariants have 
to be implemented by the SAT or projection method 
via a split form of the inviscid flux derivatives 14, 15 (en- 
tropy splitting 10 ) . This will be presented in Miiller 
and Yee 20 . 

The application of the spatial discretization of the 
perturbed Euler equations in transformed coordinates 
(15) results in a semi-discrete system of nonlinear 
ODEs 

f - R(U), (18) 

where U is the vector of the difference approxima- 
tions U' k and R is the vector of two-dimensional 

spatial difference operators operating on and 

with each element = — Jy* 4* ^ 

The periodic boundary conditions in the circumfer- 
ential n-direction in our application allow the use of 

/ g\ 

the standard sixth-order difference operator . The 
ODE system (18) is solved by the 4th-order classical 
Runge-Kutta method with a CFL condition of 1.783. 
To minimize the loss of the global order of accuracy, 
we impose the time-dependent physical boundary con- 
ditions (the in-going Riemann invariants in our case) 
only after the completion of the full step of the Runge- 
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Kutta method. The Q boundary difference operator is, 
however, employed at every stage of the Runge-Kutta 
method. Other remedies are discussed in 17 19 . 

For long time wave propagation of nonlinear sys- 
tems, even with the absence of shock waves and/or 
steep gradients, spurious high frequency oscilla- 
tions are generated by non-dissipative central spatial 
schemes. To suppress these spurious oscillations, a 
modified version of the Yee et al. 9 11 high order ar- 
tificial compressibility method (ACM) filter scheme 
is used. The original ACM filter scheme involves 
high order non-dissipative spatial base schemes and 
a c4aracteristic-based ACM filter. The ACM filter 
consists of a product of numerical dissipations and an 
ACM sensor. The ACM sensor consists of a prod- 
uct of a sensor coefficient and a gradient-like detector. 
Here, the limiter of the numerical dissipation, and the 
Roe’s averaged states are not used. For low Mach 
number CAA, the limiter which is designed for cap- 
turing discontinuities might not be necessary, unless 
shear layers, high gradient regions and/or boundary 
layers are present. 

At the completion of a full step of the Runge-Kutta 
method, the numerical solution U'"^ 1 is filtered by 

U'y 1 = U^-A tjj.k [D £ U' + D„ U']"* 1 . (19) 


k represents the fourth-order numerical dissipa- 
tion in the ^-direction defined by 

D k v' jk = K6 £ R £ C £ |A £ |5|R £ - 1 6 £ U', fc 


where 6$ is the difference operator given by <5^U' k = 
. Applying <5^ four times yields the 


U J 


j + I/2,fc- 


_T y 


standard central difference approximation - k — 


U 


j+2.k “ ^U' + 1 + ^U' j. — 4U'_ 1)A; -r 

2^(x(j&£),y(kAr)),t). Here UA involves 
differences of average states. The columns of are 
the right eigenvectors of the Jacobian matrix and 
may be found in Vinokur and Yee 11 . The eigenval- 
ues of define the diagonal matrix = diag(u^ — 
c^, ii£ + c^) , where u $ = uJ~ l £ x + vJ~ l £ y , and 

= Cs /( + (J~ l £i/) 2 • Instead of the Roe- 
average, the R,c and are evaluated at the arithmetic 
average. Here the ACM sensor is the product of a 
sensor coefficient k and a gradient-like detector. A 
simplified form of the the gradient-like detector is ex- 
pressed by the diagonal matrix Cc. The Ith diagonal 
element c l of is defined by 


+ u' 


of 


where 9, 


c j + l/2,Jk “ max (^j,fc ? ^j + l.Ar) , 
1 \ a 'j + l/-2 iJ 1 


(20) 


To avoid 


J ' k max (i Q i + i/a,fcl + l Q Ul/^ J' 10_7 )‘ 

an extra logical statement, the quantity 10~ 7 is used 


to avoid division by zero. a 


i 

>+1/2 ,* 


is the wave 


strength of the Ith wave, i.e. the Ith component of 


u ;.*- f . k is defined analogously. 

For our numerical experiments, selected nonzero pos- 
itive ACM sensor coefficient k less than or equal to 
0.10 are evaluated. 

Although Muller and Yee 20 performed numerical 
studies on a similar almost circular Kirchhoff vortex 
model, aside from application specific, the present 
simulation differs from Muller and Yee in two ways. 
The fourth-order classical Runge-Kutta temporal dis- 
cretization and the full ACM sensor (k and C^) are 
evaluated. 


Numerical Results 

We consider an elliptic generator casing with radius 
R = 2m, e = 0.00125, ft = 82.5V The stagnation 
conditions are po = 13^, cq = 330 y. Thus, the 
Helmholtz number becomes H — kR — 2QR/cq — 1 . 
The time step is chosen as Af = 0.15. A polar grid 
mildly stretched near the generator casing in the radial 
direction, and uniform in the circumferential direction 
of 129 x 24 is used. As presented later, a grid re- 
finement study indicates that this grid is sufficient to 
capture the flow physics accurately. 

We consider the elliptic generator casing started in- 
stantaneously. The initial conditions are p f = u f = 
v r — p f = 0, except for the circle r — R, where the 
analytical solution is prescribed. At t = 30, the gen- 
erator has rotated 7.5 radians. Figures 3 and 4 show 
the numerical results without and with characteristic- 
based filter, respectively. Without filter, the numerical 
solution is polluted by spurious oscillations. These are 
eliminated by the filter. The wave front started at 
r = R at t — 0 has reached r ^ 16 at t = 30. The 
numerical solution with filter is compared with the an- 
alytical solution along the positive x-axis at t = 30 in 
Fig. 5. Between x — 1 and x « 14, the analytical 
and numerical solutions are in close agreement. At 
x « 16, we see that in general the wavefront cannot 
match the infinitely long rotating generator solution, 
because the instantaneously started generator has zero 
acoustic pressure downstream of the wavefront. The 
discrepancies, therefore, have physical reasons. 

At t — 75, the generator has rotated 18.75 radians , 
Now*, the wavefront is at r « 38.5. If the wavefront 
were reflected at the farfield r = 32.1875, we should 
see it at r w 25.875. The computed acoustic pressure 
along the positive x-axis shown in Fig. 6 is in ex- 
cellent agreement with the analytical solution, if the 
characteristic-based ACM filter with a uniform sen- 
sor coefficient k = 0.05 is used. Without the filter, 
spurious oscillations occur. The w r avefront generated 
at t = 0 by the elliptic generator casing has passed 
through the farfield without visible reflection. The 
quadrupole structure of the acoustic pressure is cor- 
rectly recovered by the high order Q operator with 
the characteristic-based ACM filter of k — 0.05. See 
Figs. 7 and 8 for the comparison. 
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Aooutdc prtuura, H ■ k R « 1. «-0 00125 roi»t»d angl* ■ 7 5 r*d 



Fig. 3 Computed acoustic pressure contours with- 
out filter for instantaneously started generator 
with H = 1, € = 0.00125, rotated angle = 7.5rad. 


Fig. 5 Comparison of acoustic pressure for instan- 
taneously started generator with H = 1, e = 0.00125, 
rotated angle = 7.5rad. 



Fig. 4 Computed acoustic pressure contours with 
filter (* = 0.05) for instantaneously started genera- 
tor with H = 1, e — 0.00125, rotated angle = 7.5rad. 

We note that the multistage method used in Miiller 
and Yee 16 for time integration yields almost identical 
results as the classical Runge-Kutta method employed 
here. Since both ODE solvers coincide for linear prob- 
lems and we are solving essentially a linear problem, it 
is not surprising that the results with the two different 
time integrations differ 3 to 4 orders of magnitude less 
from each other than from the analytical solution. 

A grid refinement study was performed with 65 x 15, 
129 x 24 and 257 x 42 grids keeping the CFL num- 
ber constant. The maximum CFL numbers in the 
radial and circumferential directions were about 1 and 
0.2, respectively. The high order Q operator with 
the characteristic-based ACM filter of «: = 0.05 was 
used. The computed acoustic pressure along the pos- 
itive z-axis after time t = 75 in Fig. 9 indicates grid 
convergence. The ^-norm of the error of p ' along the 
positive z-axis, i.e. the deviation of the computed from 


Fig. 6 Comparison of acoustic pressure for instan- 
taneously started generator with 7i = 1 , e = 0.00125, 
rotated angle = 18.75rad. 

the analytical acoustic pressure there, has a conver- 
gence rate of 4.7 from the course to the medium grid 
and of 1.9 from the medium to the fine grid. The error 
of p ' along the positive z-axis is shown in Fig. 10. The 
reduction in convergence rate seems to be due to the 
boundary treatment, because the error is largest near 
the boundaries. A similar drop of convergence rate 
is observed for the conservative variables in the whole 
domain as well. 

The influence of the gradient-like detector of the 
ACM sensor and its sensor coefficient k is investigated. 
We consider the 129 x 24 grid for the same test case 
as above. This test case consists of smooth flow and it 
is not a very long time wave propagation simulation, 
consequently, only a small sensor coefficient is needed. 
With a small sensor coefficient, the role of the gradi- 
ent detector is diminishing. This is due to the fact 
that the gradient-like detector is designed to minimize 
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Acoutfjc pTwiur*. Hikfl. 1 1-000125. ro*at*d - 18 T5 'Mi 


Acousflc |X*s«ur« H « k R - 1 «-0 00125. rotated »ngt# ■ 18 75 r*t 



Fig. 7 Acoustic pressure contours with filter 
(k — 0.05) for instantaneously started generator 
with Ti = 1, € = 0.00125, rotated angle = 18.75rad. 


Acoustic p»**iur» H ■ k R * t. 00125 r o1»1»d angle « 18 75 rad 



Fig. 8 Exact acoustic pressure contours for gen- 
erator with Ti = 1, e = 0.00125, rotated angle 
= 18.75rad. 

numerical dissipation away from shocks, shears, high 
gradient and steep boundary- layer flows with the use 
of a larger sensor coefficient. However, for very long 
time wave propagation problems and/or with unknown 
solution behavior, the full ACM sensor or wavelet sen- 
sor (cf. Sjogreen and Yee 21 ) is necessary, even for 
low Mach number CAA (cf. Yee et al. 10 ). As ex- 
pected, the difference in results with and without the 
gradient-like detector (C$ of (20) vs. = I) are not 
pronounced for the very small sensor coefficients stud- 
ied. For k = 0.025, the error of the acoustic pressure 
is rougher with the gradient detector than without, 
cf. Fig. 11. For k = 0.05, the error of the acous- 
tic pressure with the gradient detector and without 
behave similarly. However, for k — 0.1, the error of 
the acoustic pressure is larger without the gradient 
detector than with the detector, cf. Fig. 12. In this 
case, the gradient-like detector is performing its duty 


x 10* 



Fig. 9 Grid refinement study of acoustic pressure 
for instantaneously started generator with H = 1, 
e = 0.00125, rotated angle = 18.75rad. 


Em* in acoustic prossuo. H -k fl ■ 1, r-0 00125. rotated angle - 18.75 red 



Fig. 10 Grid refinement study of error in acoustic 
pressure for instantaneously started generator with 
7Y = 1, 6 = 0.00125, rotated angle = 18.75racf. 

to minimize the use of numerical dissipation where it 
is not needed (in this case, away from spurious oscil- 
lations). Thus, the results seem to be less sensitive 
to the filter coefficient k , when the gradient detector 
switch is used. 

The result of Yee et al. 10 indicated that the use 
of entropy splitting of the inviscid flux derivative can 
minimize the use of numerical dissipation. For the 
present model problem, a similar conclusion as Yee et 
al. can be drawn. With entropy splitting, w r e can ob- 
tain reasonable results even without filter. Figures 13 
and 14 compare the results without and with entropy 
splitting for p* along the positive z-axis. The improve- 
ment with entropy splitting is clearly visible. More 
results and the formulation of the entropy splitting 
for low Mach number aeroacoustics will be reported in 
Muller and Yee 20 . 
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Fig. 11 Influence of the ACM sensor with n = 0.025 
on error of acoustic pressure for instantaneously 
started generator with Ti = 1, e = 0.00125, rotated 
angle = 18.75rad. 


E W in *C 0 urtC priuur* «>0 001 75. rol*t*d *ngt« - 1*75 f*t 



Fig. 12 Influence of the ACM sensor with k = 0.1 
on error of acoustic pressure for instantaneously 
started generator with Ti — 1, € = 0.00125, rotated 
angle = 18.75rad. 

Conclusions 

The periodic deformation of electrical generators by 
the rotors’ magnetic fields is modelled as a rotating 
ellipse. The normal velocity of the elliptic genera- 
tor surface provides the boundary condition for the 
acoustic field. For an almost circular generator sur- 
face, the Helmholtz equation governing the acoustic 
pressure is solved analytically in terms of the Hankel 
function of second order. Thus, the sound emission of 
rotor induced deformations of generator casings can be 
described analytically. The resulting model problem 
provides a good test case for the evaluation of CAA 
numerical algorithms. 

Although the numerical simulation of the resulting 
almost circular Kirchhoff vortex model is similar to the 
study in Muller and Yee , 16 a different form of the filter 


AeouCtC pm«n. H»kRi 1, «-<J 00125 ro«**d M"V* . 18 75 r*d 



Fig. 13 Influence of no entropy splitting on acous- 
tic pressure for instantaneously started generator 
with H — 1, c = 0.00125, rotated angle = 18.75rad. 


Acoutfc prMKjr*. H - K R - l.c-0.00125. . 1* 75 -*d 



Fig. 14 Influence of entropy splitting on acoustic 
pressure for instantaneously started generator with 
H — 1, (. — 0.00125, rotated angle = 18.75rad. 

and temporal discretization are evaluated, and grid re- 
finement study is also included. Results reveal a very 
accurate simulation with a fairly coarse grid for yet an- 
other important CAA physical application. For CAA 
flow problems of this type, especially for long time 
wave propagations, appropriate filter applied to the 
entropy splitting perturbation form of the Euler equa- 
tions is essential to achieve stability and high accuracy 
of high order non-dissipative spatial base schemes. 
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